Nuclear masses set bounds on quantum chaos 
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It has been suggested that chaotic motion inside the nucleus may significantly limit the accuracy 
with which nuclear masses can be calculated. Using a power spectrum analysis we show that the 
inclusion of additional physical contributions in mass calculations, through many-body interactions 
or local information, removes the chaotic signal in the discrepancies between calculated and measured 
masses. Furthermore, a systematic application of global mass formulas and of a set of relationships 
among neighboring nuclei to more than 2000 nuclear masses allows to set an unambiguous upper 
bound for the average errors in calculated masses which turn out to be almost an order of magnitude 
smaller than estimated chaotic components. 
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The importance of an accurate knowledge of nuclear 
masses to understand the processes occurring in astro- 
physical phenomena has been abundantly stressed 0. 
Though great progress has been made in the challeng- 
ing task of measuring the mass of exotic nuclei, theo- 
retical models are still necessary to predict their mass in 
regions far from stability |2| ■ Advances in the calculation 
of atomic masses have been hampered by the absence of 
an exact theory of the nuclear interaction and by the dif- 
ficulties inherent to quantum many-body calculations, so 
diverse models which attempt to bring forth the funda- 
mental physics of the atomic nucleus have been devised. 
The simplest approach is that of the liquid-drop model 
(LDM) 3]. It incorporates the essential macroscopic 
terms, which means that the nucleus is pictured as a very 
dense, charged liquid drop. Including the discrete char- 
acter of the nucleons and their basic interactions requires 
more sophisticated treatments. The finite-range droplet 
model (FRDM) 0, which combines the macroscopic ef- 
fects with microscopic shell and pairing corrections, has 
become the de facto standard for mass formulas. A mi- 
croscopically inspired model has been introduced by Du- 
flo and Zuker (DZ) 0|(|0 with positive results. Finally, 
among the mean-field methods it is also worth mention- 
ing the Skyrme-Hartree-Fock approach 00. These mass 
formulas can calculate and predict the masses (and often 
other properties) of as many as 8979 nuclides 0, but 
it is in general difficult to match theory and experiment 
(for all known nuclei) to an average precision better than 
about 0.5 MeV 0. This minute quantity, corresponding 
to less than a part in 10 5 of the mass of a typical nucleus, 
still represents a significant fraction of the energy released 
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in nuclear decays, strongly affecting the extrapolations of 
proton and neutron separation energies required in astro- 
physical processes . 

Can these deficiencies of the mass formulas have a 
chaotic origin? It was recently suggested that there might 
be an inherent limit to the accuracy with which nuclear 
masses can be calculated 0, due to the presence of 
chaotic motion inside the atomic nucleus . The nature 
of quantum chaos is still an open question and diverse 
points of view coexist in the literature [l2L H^ |. Clas- 
sical chaotic behavior has been understood in terms of 
deterministic equations whose time evolution, however, 
has a sensitive dependence on their initial conditions. A 
different kind of unpredictability arises from the quan- 
tum nature of atoms, a behavior codified in Heisenberg's 
uncertainty principle. Yet, when these two sources of un- 
predictable features combine, the result is very difficult to 
fathom. It would appear that quantum uncertainty can 
smooth the characteristic signals of chaos, but experi- 
ments with Rydberg atoms in strong magnetic fields [l4| 
and quantum dots [15j clearly display its footprints. 

The presence of chaotic motion in highly excited nu- 
clear systems can be gauged through the statistics of 
the high-lying energy levels or resonances [l6l IrH Il9| 
and their comparison with Random Matrix Theory 
(RMT) (2jj. It has been convincingly shown that the 
power spectrum (Fourier transform squared) of the fluc- 
tuations of a chaotic quantum energy spectra is charac- 
terized by 1// noise 18]. The basic corollary of RMT [3 
is that near or around the neutron absorption threshold, 
predictability is hopeless and only a statistical analysis 
makes sense. The possible existence of a remnant of chaos 
at the level of the ground state 0| was addressed by 
Bohigas and Leboeuf ^lj from a novel perspective. The 
errors among experimental and calculated masses in 
are interpreted in terms of two types of contributions. 
The first one is associated with a regular part, related to 
the underlying collective dynamics (LDM), plus a shell- 
energy correction, while the other is assumed to arise 
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from some inherent dynamics, possibly higher-order in- 
teractions among nucleons [ll|, that lead to chaotic be- 
havior. 

In this context it is worth to analyze the differences be- 
tween measured and calculated masses using well-known 
mass models, looking for signals of quantum chaos and, 
if these are found, to figure out if such chaos is inherent 
to nuclear masses or an artifact of the approximations in 
the many-body theories. It is our purpose to also test 
the assertion that a lower limit on how accurately nu- 
clear masses can be calculated may have already been 
reached [ljj. To estimate the best available accuracy, 
we use a different approach to calculate known nuclear 
masses. Besides the "global" formulas of which FRDM 
has become the standard, there are a number of "local" 
mass formulas. These local methods are usually effec- 
tive when we require the calculation of the masses of 
a set of nuclei, which are fairly close to other nuclei of 
known mass, exploiting the relative smoothness of the 
masses M(N, Z) as a function of neutron (N) and proton 
(Z) numbers to deduce systematic trends. Among these 
methods there are the equations connecting the masses 
of particular neighboring nuclei known as the Garvey- 
Kelson (GK) relations pH 12*2^. These relations do not 
involve free parameters and can be derived from a sim- 
ple nuclear-model picture. Strictly speaking, they do not 
yield an independent calculational tool, but they do pro- 
vide strong indication that a large fraction of the mass 
values have a smooth and regular behavior. In this in- 
terpretation, the GK relations can be viewed as a simple 
methodology to estimate nuclear masses from those of its 
neighbors. 

In this Letter we carry out a systematic analysis of 
mass errors using different global mass formulas, supple- 
mented with a study of the estimates arising from the 
GK relations which we use to set an upper limit on the 
intrinsic precision with which masses can be estimated 
when enough local information is available. 

There are differen t types of GK relations [2j| , and the 
two best known are [HlHIll 

-M(N+1, Z-2)+M(N+l, Z)-M(N+2, Z-l) (1) 
+M(N+2,Z-2)-M[N, Z)+M{N,Z-1) = 0, 

M(N+2, Z)-M(N, Z-2)+M(JV+l, Z-2) (2) 
-M(N+2,Z-1)+M[N,Z-1)-M(N+1,Z) = 0. 

These simple equations are based on the independent- 
particle shell model and, furthermore, constructed 
such that neutron-neutron, neutron-proton, and proton- 
proton interactions cancel. Both GK relations provide 
an estimate for the mass of a given nucleus in terms of 
five of its neighbors. This calculation can be done in six 
different forms, as we can choose any of the six terms 
in the formula to be evaluated from the others. Using 
both formulas, we can have a maximum of 12 estimates 
for the mass of a given nucleus, if the masses of all the 
required neighboring nuclei are known. Of course, there 
are cases where only 11 evaluations are possible, and so 



TABLE I: The mass rms deviations <7 rms , in keV, for LDM, 
FRDM, DZ, and different GK calculations. 





LDM FRDM 


DZ 


GK-1 


GK-4 GK-7 GK-12 


A > 16 
A>60 


3211 653 
3177 529 


362 
320 


171 
115 


131 112 86 
100 88 80 



on. About half of all nuclei with measured masses [24J 
can be estimated in 12 different ways and, in all cases, 
our estimate corresponds to the average value. To our 
knowledge, the systematic application of GK relations in 
this extended fashion is new. 

Since from this simple method we obtain a specific 
mass value entirely determined by the behavior of the 
neighboring nuclei, we may conclude that these values 
represent a uniform and smooth component of the nu- 
clear masses. In this sense they provide an upper esti- 
mate of the intrinsic precision limit with which they can 
be calculated, when enough local information is available. 
Fluctuations on top of the regular behavior predicted by 
different mass formulas and the GK relations are ana- 
lyzed in the second part of this paper. 

In what follows we compare the mass deviations found 
in three of the global methods (LDM, FRDM, DZ) and 
those in the GK estimates made above. The correspond- 
ing root-mean-square (rms) deviations for N nuclei 
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are displayed in Table [IJ which also shows deviations for 
the smaller samples GK-rt which involve the application 
of at least n GK relations. Note the decrease in the er- 
rors from LDM to DZ, as a consequence of the inclusion 
of shell corrections (FRDM) and terms that mimic two- 
body interactions (DZ). A similar effect is also seen in the 
GK estimates with an increase in precision, proportional 
to the number of GK relations applied. We conclude from 
this analysis that, if an intrinsic precision limit exists at 
all in the calculation of nuclear masses, it is smaller than 
100 keV. We remark that the uncertainty range implied 
by these calculations is consistent with recent estimates 
of the effect of statistical fluctuations of high- lying config- 
urations near the neutron threshold on the ground-state 
energy [25j. We also stress that experimental errors are 
lower than the estimates shown in Table I. For example, 
the average experimental error is cr CX pt ~ 20 keV for the 
nuclei in GK-12. 

Another relevant question is whether nuclear masses 
far from stability can be predicted with a similar ac- 
curacy. In this respect, the GK relations cannot by 
themselves offer a simple answer, but large-scale shell 
model calculations shed some light on this important is- 
sue [2(| . Sophisticated new calculations indicate that 
predictive power seems robust against long-distance ex- 
trapolations. Masses of 67 light nuclei in the fp shell, 
many of them unstable, were calculated with an average 



error of 215 keV |27|. Errors come mostly from isospin 
violation, and do not increase far from stability. Large- 
scale shell model calculations for the N = 126 Po-Pu 
isotones [2^ find errors in the binding energies smaller 
than 50 keV along the whole chain, with no increase for 
unstable nuclei, and "imply a high predictive power for 
ground-state binding energies beyond the experimentally 
known nuclei." 

We now proceed to consider the presence of quantum 
chaos and analyze to this end the statistical characteris- 
tics of the mass fluctuations present in the different mod- 
els. A direct 2D study in the N-Z plane exhibits clear evi- 
dence of the correlations between mass errors in neighbor- 
ing nuclei [2|j which, proceeding through the four models 
(LDM, FRDM, DZ, and GK), decrease in size. It is, how- 
ever, difficult to perform a systematic analysis of them 
because of the irregular form of the nuclear-data chart, 
and because cuts along fixed TV, Z, or A lines have a small 
number of nuclei, making it difficult to extract unequiv- 
ocal conclusions [29|,|3(j. To overcome these difficulties, 
we organize all nuclei with measured masses by ordering 
them in a boustrophedon single, ID list |3l| . numbered as 
follows. Even- A nuclei are ordered by increasing N — Z, 
while odd- A ones follow a decreasing value of N — Z. Fig- 
ure 1 displays the mass differences plotted against the or- 
der number, for the four sets of calculations. The LDM 
errors are conspicuously large around shell closures, as 
expected. The presence of strong residual correlations in 
these deviations is also apparent from Fig. 1. Regions 
with large positive or negative errors can be clearly seen 
in FRDM. The distribution of the DZ errors is closer 
to the horizontal axis, and the correlations are less pro- 
nounced, although not completely absent. Finally, the 
remaining GK discrepancies are small and highly fluctu- 
ating, without any remarkable feature except for larger 
errors for A < 60. 

Since the ordering provides single- valued functions, the 
discrete Fourier transforms 

1 Af^P-Mf f-2mjk\ fA . 

Ffc = ^?^T^ exp (^v-)- (4) 

can be evaluated. The corresponding spectral distri- 
butions can then be fitted to a power law of the form 
|-F(a;)| 2 ~ uj" 1 with uo — k/N. The squared amplitudes 
\Fk\ 2 are shown in Fig. 2 where the straight lines are fit- 
ted slopes, which in a log-log plot correspond to the value 
of to in the corresponding power law. While fluctuations 
are large, the results are remarkable. The gradual van- 
ishing of the slope in the different calculations from top 
to bottom is evident from the figure. 

The fluctuation study allows the application of results 
well known from noise analysis. Signals with strong cor- 
relations over a long range tend to have a Brownian mo- 
tion behavior, where their power spectrum approach the 
power law |-F(a-i)| 2 ~ oj m , with m = —2, where u> is the 
frequency of the Fourier component. At the other end 
of the correlation spectrum is white noise, where there is 
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FIG. 1: Differences, in MeV, between the measured 
masses |24j and those obtained in LDM, FRDM, DZ, and 
GK, plotted against the order number j. 



no remaining coherence, which implies a flat frequency 
dependence m = 0, that is, all frequencies have equal 
weight. Intermediate values of the slope m correspond 
to transitional regimes. The first three sets of mass de- 
viations in Fig. 2 give rise to power laws with 777 ranging 
from —1.8 to —0.5, consistent with a slow transition from 
substantial correlations to smaller (but still significant) 
ones for the associated mass fluctuations. The chaotic- 
ity discussed in Ref.JO, according to the m =-1 criteria 
put forward in Ref . [l8| , seems indeed to be present in the 
FRDM deviations, while it tends to diminish in the mi- 
croscopic DZ calculations. In contrast, we find that the 
GK analysis is essentially consistent with white noise. 

In summary, a careful use of several global mass for- 
mulas and a systematic application of the Garvey-Kelson 
relations imply that there is no evidence that nuclear 
masses cannot be calculated with an average accuracy 
of better than 100 keV. While mass errors in mean-field 
calculations like the FRDM behave like quantum chaos, 
with a slope in their power spectrum close to —1, mi- 
croscopic models' results correspond to smaller slopes. 
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Finally, for the local GK relations the remaining mass de- 
viations behave very much like white noise. These results 
seem to confirm that the chaotic behavior in the fluctua- 
tions arises from neglected many-body effects. The more 
physical information is included in the model, through 
many-body terms or by means of local information which 
measure the regular and smooth components in the mass 
systematics, the less presence of chaos is observed. The 
local calculations could in principle be reproducing the 
chaotic dynamics, which would imply that no a priori 
limit exists on predictability. A different interpretation 
is that the errors (that simulate chaotic behavior) are 
simply removed by the many-body effects. In any case, 
our results strongly suggest an upper bound of approxi- 
mately 100 keV as precision limit for the calculation of 
nuclear masses. Although the presence of (chaos-related) 
unpredictability at such small scale cannot be ruled out 
by our study, we believe this is encouraging news in the 
quest for reliable predictions of nuclear masses. 



FIG. 2: Squared amplitudes of the Fourier transforms of the 
mass differences obtained in LDM, FRDM, DZ, and GK, plot- 
ted against the frequency lu = k/N. 
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